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ABSTRACT 

It appears that the Lyman a forest is becoming thick at a redshift of about 5.5, cutting off the higher 
redshift intergalactic medium from view in neutral hydrogen. However, the effects of star formation at 
higher redshift are still readable in the intergalactic metal lines. In this paper I use observations of 32 
quasars with emission redshifts in the range 2.31 to 5.86 to study the evolution of the intergalactic metal 
density from z = 1.5 to z = 5.5. The C IV column density distribution function is consistent with being 
invariant throughout this redshift range. From direct integration, I determine ficiv to be in the range 
(2.5 — 7) x 10~ 8 and Qsiiv in the range (0.9 — 3) x 1CV 8 between z — 1.5 and z = 5. The metallicity at 
z = 5 exceeds 3.5 x 10 -4 , which in turn implies that this fraction of the universal massive star formation 
took place beyond this redshift. This is sufficient to have ionized the intergalactic medium. 

Subject headings: early universe — intergalactic medium — quasars: absorption lines — galaxies: 
formation 



1. INTRODUCTION 

The epoch of reionization is the ultimate depth to which 
we can directly investigate neutral hydrogen in the in- 
tergalactic medium (IGM). Whether or not it signals the 
onset of reionization, the increasing optical depth of the 
Lyman a forest at redshifts above 5 (Becker et al. 2001; 
Djorgowski et al. 2001; Fan et al. 2000b, 2001c; Songaila 
et al. 1999) severely limits our ability to study the IGM 
with neutral hydrogen at these times. However, we can 
continue to study the metal lines in the quasar spectra, 
and in particular the C IV 'forest'. These absorption lines 
can be used to infer the history of early star formation by 
mapping the evolution of the universal density of metals 
with redshift. The contribution to the universal metal- 
licity arising from the diffuse IGM has been extensively 
analyzed at redshift z ~ 3 (Cowie et al. 1995, Tytler et al. 
1995, Songaila 1997, Songaila 1998, Ellison et al. 2000) 
and some information is also available from investigations 
of O VI at lower redshift (Buries & Tytler 1996, Tripp, 
Savage & Jenkins 2000). However, to date there has been 
no systematic analysis of the variation of the intergalactic 
metal density as a function of redshift. 

Most interpretations of the origin of the z = 3 metal- 
licity have focussed on ejection or stripping of material 
from small galaxies that formed at higher redshifts (z > 5) 
(Gnedin & Ostriker 1997, ; Gnedin 1998; Madau, Ferrara 
& Rees 2001). It is also possible that the metals may 
originate in genuine Population III star formation at very 
high redshift (e.g., Carr, Bond & Arnett 1984; Ostriker 
k Gnedin 1996; Haiman & Loeb 1997; Abel et al. 1998). 
In these scenarios most of the IGM metals would be in 
place at the highest currently observable redshifts. More 
recently, Steidel has suggested that the observed metals 
may instead originate in extremely fast winds from Ly- 
man break galaxies at z — 3 (e.g. Pettini 2001). Such a 
model may already encounter problems in explaining the 
kinematic quiescence of the C IV lines (Rauch, Sargent 
& Barlow 2001) but it can also be tested by determining 
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whether the metallicity changes with redshift. 

This paper analyses the distribution of C IV and Si IV 
ions in the IGM over a wide range of redshifts from z = 1.5 
to z — 5.5 using newly-discovered Sloan Digital Sky Sur- 
vey (SDSS) quasars (Anderson et al. 2001; Fan et al. 1999, 
2000a,b, 2001a,b; Zheng et al. 2000) to provide a high red- 
shift sample. I use the methodology previously applied 
at z ~ 3 (Songaila 1997) of directly integrating the ob- 
served ion column densities in Lya forest clouds to obtain 
strict lower limits on f2i on , which, with plausible assump- 
tions on ionization corrections and abundance patterns, 
can provide a strict lower limit on r2 me tais at a number of 
redshifts. 

New observations of z = 2 — 3 quasars with the HIRES 
spectrograph (Vogt 1994) on the Keckl telescope and of 
z — 4 — 5.8 quasars with the ESI instrument (Sheinis et 
al. 2000) on Keckll are described in section 2. In section 3 
I construct column density distribution functions for C IV 
and Si IV from z = 1.5 to z = 5.5 and perform the inte- 
grations. Conclusions are given in section 4. 

2. OBSERVATIONS 

The data used in the paper comprise HIRES observa- 
tions of 13 quasars reported in Songaila (1998) together 
with HIRES and ESI observations of a further 19 quasars. 
The reduction of the HIRES (lower redshift) observations 
follows the procedures given in Songaila (1998). The 
fainter, high redshift, quasars were observed with the ESI 
instrument in echellette mode. (See Songaila (2001) for a 
list of quasars observed with ESI.) The resolution in this 
configuration is comparatively low, ~ 5300 for the 0.75" 
slit width used, but the red sensitivity is high and the 
wavelength coverage is complete from 4000 A to 10, 000 A. 
At the red wavelengths, precision sky subtraction is re- 
quired since the sky lines can be more than two orders 
of magnitude brighter than the quasars. In order to deal 
with this issue, individual half-hour exposures were made, 
stepped along the slit, and the median was used to pro- 
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vide the primary sky subtraction. The frames were then 
registered, filtered to remove cosmic rays and artifacts, 
and then added. At this point a second linear fit to the 
slit profile in the vicinity of the quasar was used to re- 
move any small residual sky. The observations were made 
at the parallactic angle and flux calibrated using obser- 
vations of white dwarf standards scattered through the 
night. These standards were also used to remove telluric 
absorption features in the spectra, including the various 
atmospheric bands. The final extractions of the spectra 
were made using a weighting provided by the profile of 
the white dwarf standards. Wavelength calibrations were 
obtained using third-order fits to CuAr and Xc lamp ob- 
servations at the begining and end of each night, and the 
absolute wavelength scale was then obtained by register- 
ing these wavelength solutions to the night sky lines. The 
wavelengths and redshifts are given in the vacuum helio- 
centric frame. 

I next constructed the sample of all C IV and Si IV 
aborption lines which lie more than 50 A longward of each 
quasar's Lyman alpha emission and hence clearly outside 
the Lyman alpha forest region. The doublets were found 
by inspection of the spectra, and confirmed by consistency 
of the column density and velocity structure in the two 
members of the doublet. Systems within 5000 km s _1 
of the quasar's emission redshift were excluded to avoid 
contaminating the sample with proximate systems whose 
ionization may be dominated by the quasar itself. The 
final C IV sample consists of 367 C IV doublets, with 
1.74 < z < 5.29 and 12.0 < log(A c w ) < 14.9. The 
Si IV sample contains 109 systems between z = 1.78 and 
z = 5.29, with 12.0 < log(iV si IV ) < 14.8. 
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Fig. 1. — Column density distribution of C IV absorbers for 
the present samplc(fillcd symbols) in the redshift interval 2.90 < 
z < 3.54; / is the number of systems per column density interval 
per unit redshift path. The bin size is 10 0,3 ./V cm -2 and the error 
bars are ±1<t based on the number of points in each bin. The solid 
line is a power law of the form f(N)dN = BN~ a dN with index 
a = 1.8±0.1 established by a maximum likelihood fit to all systems 
with AT(CIV) > 10 cm" . The open symbols are the column den- 
sity distribution of Ellison ct al. (2000) for C IV absorbers in the 
quasar Q1422+231 over the same redshift interval (not corrected 
for incompleteness at low N). The dot-dash line is Ellison et al.'s 
power law fit with index 1.44 ± 0.05 to TV(CIV) = 10 12 cm~ 2 . 
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Fig. 2. — Column density distribution of C IV absorbers as in 
Figure 1 for various redshift ranges. In each redshift interval, the 
solid line is the fiducial power law f(N)dN = BN- a dN with the 
normalization and index a = 1.8 established from Figure 1. 

In all cases, ion column densities were determined by 
fitting up to ten Voigt profiles to each redshift system, 
defined to be all absorption near a given fiducial redshift 
with gaps in velocity space of no more than 50 km s _1 . 
The column density at a given redshift is then the total 
column density of all such components. In all but a few 
systems, individual lines were unsaturated, so that fitted 
column densities are insensitive to &-value. For the ESI 
spectra, where the lower resolution observations make this 
issue more critical, I adopted in doubtful cases a column 
density based on assuming that the weaker member of the 
doublet was unsaturated. This provides a minimum esti- 
mate of the column density. 

3. DISTRIBUTION FUNCTIONS 

With these samples, the column density distribution 
function f(N) was determined for each ion. f(N) is de- 
fined as the number of absorbing systems per unit red- 
shift path per unit column density, where at a given red- 
shift z the redshift path, X(z) is defined as X(z) = 
| [(1 + zfl 2 - 1] for q = 0.5 (Tytler 1987). Throughout 
the paper the distributions are calculated in c.g.s. units 
and are plotted with 1 a error bars calculated from the 
Poisson errors based on the number of systems in each bin 
(Gehrels 1986). In addition, all calculations were made 
assuming a, q n — 0.5, A = cosmology, which is a good 
approximation to the currently favored A-dominated cos- 
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mologies above about 2 = 1. 

In Figure 1 I show the column density distribution in 
the redshift interval 2.90 < z < 3.54 and compare this 
with the distribution over the same redshift range from 
an extremely deep exposure of the quasar 1422+231, as 
reported in Ellison ct al. (2000). From the figure it can 
be seen that the present data become incomplete below a 
column density of approximately 10 13 cm~ 2 in this red- 
shift range. Above this value a maximum likelihood fit 
to a power law slope gives an index of a = 1.8 ± 0.1. 
This is steeper than the Ellison et al. fit of 1.44 ± 0.05 
based on the more limited sample, but in fact provides a 
good fit to those data points at column densities less than 
10 13 cm~ 2 and to the corrected counts in 1422+231 at 
less than 10 12 cm~ 2 . I adopt this fit, with a normalization 
of log / = —12.4 at a column density of 10 13 cm -2 , as a 
fiducial to compare with the distributions in other redshift 
ranges. 

The C IV and Si IV distribution functions are shown as 
a function of redshift in Figure 2. For rcdshifts less than 
4 the data is drawn from HIRES observations of relatively 
bright quasars, and significant incompleteness occurs only 
below column densities of about 10 13 cm~ 2 . At higher red- 
shifts, where most of the observations are from ESI and 
the quasars arc fainter, simulations show there is clearly 
incompleteness even above this column density. In the 
highest redshift bin, 5 — 5.5, lying at the longest wave- 
lengths where the night sky is highest, recovery of lines 
is quite sporadic and the sample is extremely incomplete. 
Over the column density range in which incompleteness is 
not a problem, the distributions are consistent with being 
invariant throughout the range of redshifts sampled. 
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Fig. 3. — Sl(CIV) as a function of redshift plotted from the data 
of Table 1 (Ho = 65 km s _1 Mpc -1 and go = 0.5). In each redshift 
interval, the filled symbols show il(CIV) computed for all systems 
with 10 12 cm~ 2 < TV(CIV) < 10 15 cm -2 , whereas the open symbols 
are f2(CIV) computed for all systems with 10 13 cm" 2 < AT(CIV) < 
10 14 cm~ 2 . Symbols arc plotted at the average redshift for each 
bin (Table 1). Error bars are 90% confidence limits computed using 
Monte Carlo simulations. The dot-dash lines arc f!(CIV) computed 
assuming Ui, h 2 = 0.023, h = 0.65 , a C IV ionization fraction of 0.5 
and metallicities Z = 0.0001 and 0.001, respectively. 



100.00 




a 0.10 r 



0.01 I i i i i I 

1 2 3 4 5 6 

REDSHIFT 



Fig. 4.— As in Fig. 3, for Si IV. 



Sic iv and fisi iv were next calculated from the line sam- 
ples in each redshifts interval according to 

n- - H ° n 1 ' 2 ^ Nion m . a) 

"ion — " m . v " l ion 

CPcrit A A 

where p cr it = 1.89 x 10 _29 ft. 2 g cm~ 3 is the cosmo- 

logical closure density, mi on is the ion's mass, Ho = 
lOO/i km s _1 Mpc -1 , Q, m is the local matter density rela- 
tive to the closure density, dX = (1 + z) x l 2 dz is the value 
for a flat matter-dominated cosmology and AX is the in- 
tegral of X over the redshift interval. For A-dominated 
cosmologies, the values may be simply scaled to preferred 
values of tt m using the dependence in equation (1). The 
values of Qq iv and Qsi iv are plotted in Figures 3 and 
4 and tabulated in Table 1 for H = 65 km s _1 Mpc -1 
and qo = 0.5. The error bars are 90% confidence limits; 
they were computed by using the power-law distribution 
with the shape of the fiducial fit, and normalized to match 
the measured f2; on , to generate Monte Carlo simulations of 
line lists and then measuring f2; on from these simulations. 

4. DISCUSSION 

If the ionization fraction can be estimated, then f2j on 
can be translated to f2 e i ement for each element in question. 
I shall assume here that C IV and Si IV are the dominant 
ionization stages with ionization fractions of 0.5 (Songaila 
1997). This assumption is conservative in that it will, if 
anything, underestimate the resulting mctallicity since nei- 
ther C IV nor Si IV can be a much larger fraction than this; 
even in the extreme case, the metallicity will be overesti- 
mated by no more than a factor of two. The measurements 
of f^oiemcnt can then be compared with measurements of 
fib to obtain the universal metallicity contribution of the 
gas. (The inverse procedure of guessing a metallicity and 
then using this to estimate has been followed by Tripp, 
Savage & Jenkins (2000) in analyzing O VI measurements 
at low redshift. However, given that we know Qb from 
other sources, it seems more reasonable in the present case 
to use that value to constrain the metallicity.) A value of 
fib h 2 = 0.023+0.003 is adopted from measurements of the 
microwave background (Netterfield et al. 2001), which is 
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consistent with the most recent deuterium estimates (e.g., 
Tytler et al. 2001). The metallicity is quoted relative to 
solar values of 3.3 x 10 -4 for carbon and 3.3 x 10~ 5 for 
silicon (Anders & Grevesse 1989). The results of this pro- 
cedure are illustrated with the dashed lines in Figures 3 
and 4. For h = 0.65, I conclude that the carbon metallic- 
ity is 5 x 10~ 4 at lower redshifts and still exceeds 2 x 10 -4 
at z = 5. The silicon metallicity is about 2 times higher, 
3.5 x 10~ 4 at z = 5, with the assumed Si IV ionization 
fraction of 0.5. This presumably reflects the faster gener- 
ation of the alpha process elements, but the exact amount 
of the excess of silicon over carbon will of course depend 
on the relative Si IV and C IV ionization fractions. 

The minimum nature of the estimates should be em- 
phaziscd: in any given region, only observed metals are 
counted. While we believe that the bulk of the baryons are 
to be found in these clouds at these redshifts (Rauch et al. 
1997; Kim et al. 1997), the bulk of the metals may not, and 
may instead be primarily retained in the galaxies that gen- 
erated them. Furthermore, the assumption that C IV and 
Si IV are dominant ionization stages depends critically on 
the shape and intensity of the metagalactic ionizing radia- 
tion field and on the densities of the clouds. There may be 
cloud column densities and epochs for which this assump- 
tion is no longer true and where the observed ion densities 
would imply higher metallicities. Finally, the high noise 
levels at the longest wavelengths mean we are substan- 
tially undercounting the lines at the highest redshifts, and 
the turndown in the ion densities at these redshifts may 
be due to this effect rather than to a real physical change. 
Simulations of the incompleteness above z = 4.5 suggest 
that C IV line identification is substantially complete to 



log A = 13.25 at z = 4.5 and to log TV = 13.7 above z = 5, 
resulting in corrections to f2(C IV) of factors of 1.35 and 
1.67, respectively, in these two bins. 

Despite these caveats we can still draw some substantial 
conclusions. First, since total metal densities are domi- 
nated by oxygen, which we expect to scale most closely 
with silicon, we would conclude that a minimum univeral 
metallicity of about 3.5 x 10 -4 of the solar metallicity is 
in place at z = 5. This means that this fraction of the 
universal massive star formation should have taken place 
before this time, and this in turn can be translated to a 
number of ionizing photons per baryon. Since a metallicity 
of 10~ 5 gives rise to one ionizing photon per baryon (e.g., 
Miralda-Escude & Rees 1997), the metallicity measured at 
z = 5 is sufficient to have preionized the IGM. 

Second, there appears to have been little change in the 
univeral ion densities since z — 5. It is possible that this 
could be caused by matched variations in the ionization 
parameter and the metallicity. However, at least above 
z ~ 3, where the metagalactic ionizing spectrum is likely 
to be galaxy-dominated (e.g. Steidel et al. 2001), C IV 
and Si IV will be the dominant ionization stages, and it 
appears more likely that the metal densities have remained 
relatively constant during this period. Although this does 
not rule out the possibility of additional contributions from 
mechanisms such as galactic winds, it does suggest that, 
to a large extent, the metals observed in the intergalactic 
medium at z = 3 have an earlier origin. 

This research was supported by the National Science 
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98480. 
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Table 1 

flciV AND f^SilV BY REDSHIFT 
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AX 
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# Lines 


^CIV 
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# Lines 


^SilV 






C IV 


C IV 


C IV 


(xlO- 8 ) 


Si IV 


Si IV 


Si IV 


(xlO- 8 ) 


1.5 - 


2.0 


3.54 


1.869 


14 


3.74 










2.0 - 


2.5 


7.22 


2.272 


59 


6.83 


3.67 


2.243 


13 


1.65 


2.5 - 


3.0 


11.08 


2.779 


88 


4.79 


6.05 


2.806 


15 


2.90 


3.0 - 


3.5 


7.46 


3.206 


67 


4.50 


5.58 


3.192 


28 


1.13 


3.5 - 


4.0 


10.49 


3.746 


55 


6.71 


2.66 


3.700 


8 


2.91 


4.0 - 


4.5 


11.93 


4.240 


67 


7.18 


9.69 


4.277 


31 


2.97 


4.5 - 


5.0 


5.36 


4.655 


16 


2.46 


4.35 


4.642 


9 


0.85 


5.0 - 


5.5 


1.25 


5.285 


1 


0.90 


1.25 


5.285 


1 


0.17 



